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Abstract 

The density-density correlation profiles obtained superimposing absorption images from 
atomic clouds freely expanding after the release of the confining optical lattice can be 
theoretically described in terms of a generalized quantum measure based on coherent- 
like states. We show that the corresponding density patterns differ in a testable way 
from those computed using standard many-body mean values, usually adopted in fitting 
experimental data. 



1 Introduction 

A standard technique used in experiments for extracting information on the behavior of 
ultracold atomic gases trapped in optical lattices^ is based on analysis of interference phe- 
nomena (see [3]-[llj and references therein). Since the measure of relevant observables inside 
the optical lattice is problematic, the usual adopted procedure consists in the release of the 
confining optical potentials, followed by the (free) expansion of the atomic cloud up to meso- 
scopic sizes. At this point, the cloud is illuminated by a laser beam and the corresponding 
absorption image collected. The absorption process is destructive; nevertheless, many pic- 
tures can be obtained by starting each time with a new system, prepared in the same initial 
state. By superimposing the various obtained pictures, information on the atom density at 
the moment of trap release can be inferred. 

The image that is obtained by superimposing all these "photographs" is usually interpreted 
as the average of the density operator over the state of the sample. An alternate theoretical 
description is however possible. Indeed, in quantum mechanics, any averaging procedure 
obtained through a measuring process corresponds to a generalized quantum measure, i. e. a 
so-called Positive Operator Valued Measure (POVM) [T3]- [T7]. In general, the choice of the 



1 For recent reviews on the study of quantum many-body effects in these systems, see [HE]- 
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POVM to be used is suggested by the experimental evidences. In the present case of ultracold 
atoms in optical lattices, this evidence comes from experiments involving a two-well trapping 
potential [IE]. The actual data show that interference effects appear in a single absorption 
image even when the system is prepared in a totally incoherent state: more precisely, the 
interference pattern seen in single pictures seems to always conform to what is expected for 
a condensed, fixed-phase state, a state in which all atoms share the same single particle wave 
function. By constructing a POVM in terms of these coherent-like states, one finds that the 
corresponding generalized quantum measurement process leads to predictions that, at least 
in line of principle, differ from those obtained through the simple average of the density 
operator [TDjEI 

The aim of the present investigation is to analyze a possible physical scenario in which 
those differences may become visible and experimentally detectable. To this end, we shall 
study the behavior of a system of cold atoms in bichromatic optical lattices, where a second, 
low-intensity laser is superimposed to the one forming the periodic potential with the aim 
of obtaining an unbalanced filling of the lattice sites [20]- [23]. By exploiting the properties 
of the density-density correlation function [21], one can show that there are experimentally 
relevant instances in which the differences in the predictions of the two above mentioned 
theoretical interpretations can be revealed. This result has been supported by a numerical 
simulation, reproducing the situation of an actual experimental setup. We are confident that 
these results will stimulate further direct analysis and tests. 

2 Cold bosonic gases in optical lattices 

We shall study the behavior of N bosons confined in one-dimensional lattice with M sites, 
each separated by a fixed distance <ijf| In a suitable approximation, i.e. for a large enough 
inter-site barriers, their dynamics can be described by a Bose-Hubbard Hamiltonian [D[2], 

H = — J ^ b\ bj + ^ €i hi + -U ^ hi (hi — 1) , i, j = 1, 2, . . . , M , (1) 

<i,j> i i 

where < i,j > means nearest neighbor, b\, hi are creation and annihilation operator for an 
atom in site i, hi = b\i>i is the number operator on site i, while J, U and e« are parameters 

2 The standard theoretical intepretation of these experimental facts explains the appearence of the density 
interference fringes as an "emergent phenomena" [7]- [TT] ; although phenomenological in character, also this 
averaging procedure can be interpreted in terms of a quantum generalized measure, albeit through a rather 
involved, unnatural POVM. 

3 In the actual experimental setups, the bosons are really confined in an three-dimensional harmonic trap 
over which a periodic potential along say the x direction is superimposed. Since the confining potentials are 
separable, the single particle wave function describing the state of the atoms factorizes in a part depending 
only on the variable x times a piece depending on the couple (y,z). The dependence on the transverse 
coordinates is irrelevant for the dynamics in the lattice, that therefore can be effectively described by a 
one-dimensional Hamiltonian. 
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that quantify the hopping, repulsion and single-site depth energy, respectively. The operators 
b\, hi obey the standard Bose commutation relations: [S*,St] = Sij. 

The total number of particles N is conserved by the dynamics generated by ([T]). Therefore, 
the Hilbert space of the system is ( iV+ ^ f_1 ) -dimensional and can be be spanned by the set 
of Fock states, \k; N), describing the situation in which the occupation number of each site 
is fixed; the M-dimensional vector k = (ki,kz, ■ ■ ■ ,kM), with = N, represents a 

possible distributions of the iV atoms in the M sites. These states are obtained by acting 
with the creation operators on the vacuum |o); explicitly, one has 

\k ] N)^^(b\^(bl^...(b{ I ) k -\0) , (2) 



where in short k\ = k±\ k%\ . . . 

A different basis in the system Hilbert space is given by the collection of fixed-phase, 
coherent-like states, 



M 



N 



where £, ip are M-dimensional vectors whose components £j,(pj, j = 1,2,..., M, represent sets 
of real parameters such that <pj G [0, 2tt], £j e [0, 1] with ^\ 6j = 1- These states describe a 
physical situation in which all N atoms are in a coherent superposition, where measures 
the probability of finding an atom in the z-th site, while pi gives the corresponding phase. 
Only relative phases are relevant, so that one can arbitrarily fix the value of one of the tp j. 

The Fock states form an orthonormal set, while the coherent ones become orthogonal only 
in the large N limit nevertheless, they form an overcomplete set of states 



(iV + M-1)! f^dp, dyjtf-i f 1 „ /■i-6-6--eM- 2 



(4) 



By expanding a coherent state over the Fock basis, one finds: 

[ah - / M k \ 

\i;&N) = J2\ hr eik 'H H&)*)$ N )> ( 5 ) 



.3=1 



where the sum runs over all possible M-vectors k, whose components fcj obey the constraint 
ki = N. The overlap between a Fock and a coherent state is then given by: 




(h;N\£,ftN) = W-j^e** . (6) 



4 Indeed, for an M-site lattice, one finds that: (£ # iV|^, AT) = ( J2iLi V^ e * (v " _¥ '* ) ) 5 usin § the 
Cauchy-Schwartz inequality, one easily sees that the modulus of sum in the bracket is always less than one, 
unless £ = f and </? = </?'; as a consequence, its AT-th power become vanishingly small as N become large. 



3 



As well known [H [2] , the Hamiltonian (CQ) describes a cross-over between a superfluid and 
insulator phases, which becomes a true quantum phase transition, with order parameter 
depending on the ratio J/U, in the limit of an infinite number M of wells. Neglecting the 
shifts 6j, for small J/U, the ground state of the system is given by a Fock state (Mott 
insulator phase). On the other hand, for large J/U, the system shows phase coherence; all 
iV particles are in the same superposition and the ground state of (00) can be approximated 
by |£; tp; iV), with definite relative phases and occupation probabilities (superfluid phase). 



3 Many-body states 

In a typical experimental setup, the system of N atoms is first cooled to very low tempera- 
tures, of the order of few tens of nanokelvin, and then trapped in the optical lattice. Since 
measures of relevant observables directly in the lattice are difficult, indirect information on 
the dynamics of the atoms are usually obtained by switching off the periodic confining poten- 
tial and letting the atom gas expand freely up to mesoscopic dimensions. Absorption images 
of the expanded sample are then collected by shining it with a probe laser; by superimposing 
the various obtained pictures, information on the atom density at the moment of trap release 
can be inferred. 

In order to theoretically describe this process of measure, it is convenient to use a second 
quantized many-body formalism. Let us first introduce the field operator ip\x), creating 
from the vacuum an atom at position x, ^(x) |0) = \x); it can be decomposed as 

oo 

ft(x)=TTw i (x)b\, (7) 

i=l 

in terms of a complete set of single-particle wave functions Wi(x) = (x\wi) = (x\b\\0). 
Although the states \wi) = b\\0), i = 1,2 . . . , M, obtained by the action of the creation 
operators b\ on the vacuum are enough to describe an atom confined in the lattice, a complete, 
infinite set is needed to properly represent its state when it moves freely in space. These 
states are orthonormal and therefore one can invert ((7|) and write 

J! = /d* «,<«)*(,), (8) 

so that from [Sj,St] = (tUi\wj) = one recovers the standard bosonic (equal-time) commu- 
tation relation, ip(x) , ip^(x') = 6(x — x'). 

Assume now that the confining potential is released at time t — and denote by U% the 
unitary operator that evolves freely in time the initial one-particle states: 

\ Wi ,t) = U t \wi) = U t b\\Q) = 6{(t)|0) , b\(t) = u t l\u\ . (9) 
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The corresponding wave function is given by Wi(x; t) = (x\wi(t)) = (x\b\(t)\0), and coincides 
with a transformed Wi(x) under a ballistic expansion. At the moment of the release of the 
lattice, the Wi(x, 0), i = 1,2, ...,M, are wave functions localized at the lattice sites xf 
they can be identified with one-dimensional Wannier functions. After a large enough free 
expansion time t, sufficient for the clouds coming from the various sites to overlap, one finds 
through (JOJ) that these functions have a common envelope, differing only by a phase 



Wi{x,t) = |w(x,t)|e^ (x -^ )2 ; (10) 

indeed, for those times, the scale over which \wi(x, t)\ varies is larger than the product Md, 
giving the original dimension of the lattice. 

Since the dynamics is free, every particle in a many-body state will evolve independently 
with U t \ therefore, the evolution up to time t of the t — coherent state |£, <p; will simply 
be given by 

N 

|o> , (ii) 



and analogously for Fock states 



\k;N,t) 



^(O)* 1 ^))*"-..^^^) • (12) 



In this picture, only states evolve in time while operators, like ip{x), remain fixed. Since 
at each instant of time t the collection {|w i (t))}^ 1 is a complete set of single particle states 
obtained from the vacuum by the action of the creation operators b\{t), ip(x) can be equiv- 
alently decomposed as ip(x) = Y^=\ w i{ x 'i^) bi(t) for all times. Using this, one can compute 
the action of the field operator ^(x) on the two class of states, obtaining: 



i>(x)\g,<p;N,t) = VN (^VTje^Wjix,^ \g,fiN-l,t) , (13) 
ip(x)\k;N,t) = s ^^k~w j (x,t)\k u ...,k j - 1, k M \ N - 1, t) . (14) 



4 Generalized quantum measures 

In order to apply the previous formalism to the theoretical interpretation of the above men- 
tioned procedure of measuring density profiles, it is useful to recall some results deduced 
from the experiment. 

When M = 2 and the atoms in the lattice just before the release of the confining potential 
are prepared in a superfluid state described by (pJJ), the picture that is obtained at time 
t after a free expansion shows a high visibility interference pattern, with fringe spacing 
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mediated by the wave vector Q = md/ht, where m is the atom mass [18J. This is to be 
expected, since in the superfluid phase essentially all iV particles occupy the same quantum 
state. The roughness and imperfection of the interference figure in a single image, beside to 
experimental errors, has to be ascribed to the fmiteness of the particle number N. Indeed, 
in many-body physics, one can assimilate ensemble averages with mean values with respect 
to macroscopically occupied many-body states, provided the number of particles involved is 
large enough. Therefore, the larger the number N of atoms the system contains, the better 
a single absorption image will model the average rigJx, t) of the density operator at point 

n(x) = i/)*(x) i>(x) , (15) 

in the state £, <p;N,t). Using (fT3l) . one can easily compute the theoretically predicted 
density profile to obtain 

n l? (x,t) = {l(p-N,t\h(x)\l^N,t) =N\Y,Vl^ 3w ^^)\ 2 

j 

= iV| w (x,t)| 2 {l + 2^ v ^cos(^(^^)-^ + ^)} , (16) 

3<l 

where, recalling ffTUl) . the relative dynamical phases take the explicit form 6ji(x,t) = 
£t[—(x — %j) 2 + (x — xi) 2 ]; since for the lattice site positions one has Xj = jd, it can 
also be rewritten as 9ji(x,t) = Q(j — l)(x — |(j + /)), thus reproducing (for M = 2) the 
observed interference pattern mediated by the wave vector Q = md/ht. 

Since taking the picture of the expanded gas is a destructive operation (the sample is lost), 
one usually repeats many times the whole measuring process, preparing the system in the 
same superfluid initial state. In the case M — 2, all subsequent absorption pictures show the 
same interference figure (modulo experimental uncertainties), with the same fringe spatial 
positions. Therefore, by superimposing all these pictures, one obtains an average image that 
is indistinguishable from any of the single shot pictures, and therefore it is again described 
by the pattern ffTUl) . 

If instead one starts at t = with a system prepared in an incoherent state, like the Fock 
states in (J2j), with no definite relative phase relations, no interference pattern in the absorp- 
tion images is expected to appear. Surprisingly, this is not the case: in any experiment so far 
performed, single shot images always show interference figures compatible with the pattern 
in (1161) . irrespectively from the state of the sample at the moment of trap release. However, 
unless the sample is prepared in a superfluid state, the absolute spatial position of fringes in 
repeated single shot images is seen to vary randomly; as a consequence, by superimposing 
many single shot images, the interference pattern indeed completely disappears. 

The usual theoretical interpretations of these experimental facts is based on the assumption 
that the superposition of single shot images reproduces operator averages. More precisely, 
the procedure of repeatedly preparing the system at £ = in a generic many body state 
described by the density matrix p, and then superimposing all the single shot pictures taken 
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after a free evolution of the system up to time t, p i— > p(t), is commonly believed to reproduce 
the quantum average 

n p (x,t) = Tr[n(x)p(t)] . (17) 

For instance, for generic M, starting from a Fock state ([2]), p = \k;N}(k;N , in this way 
one would obtain 

M 

n % (x,T) = (k;N,t\h(x)\k;N,t) = kj\wj(x, t) \ 2 , (18) 

which shows no interference figures. 

An alternate, different theoretical formalization of the experimentally used measurement 
technique is nevertheless possible. It takes into account the previously mentioned result 
that, for M — 2, single shot images always present interference patterns compatible with 
the profile in (|16|) . irrespectively from the initial state p, and extends it to the case M > 2. 
In other terms, the procedure of "taking a photograph" of the expanded sample seems to 
select a coherent-like state |£, <p; N, t) for the state of the system, with definite amplitudes 
£i and phases <fi, and as a consequence produce the profile n^(x,t) in ( Tl6i) for the average 
densitylf] As seen in experiments, amplitudes and phases nevertheless vary from shot to 
shotE The occurrence of given values £ and ip for such parameters in a single shot will be 
determined by the initial state p. More precisely, the distribution of & and ipi over many 
absorption pictures will be determined by the probability (£, (p; N, t\p(t) £, (p; N, £), which 
gives the weight of the configuration |£, <p; N, t)(£> <p; N, t| in the expansion of p(t) in the 
basis of coherent states. 

As a result, in this scheme, the mean value of the density n(x), as for any other observable, 
is given by the sum over all possible configurations {£, 0} of the average ngJx, t) weighted 
with the above mentioned probability; explicitly, one then should write: 



n, 



(x,t) = / dfi(<p) / dp(0 (l<p;N,t\p(t)\l&N,t) n t Jx,t) , (19) 



where for the the properly normalized volume and integration measure, we have used the 
shorthand notation [cf. (|3J)]: 



N\ Jo 27T 27T Jo Jo 



(20) 

The expression ( fl9|) for the averaged density corresponds to a quantum mechanical gen- 
eralized measure. It can be described by an operation of trace of h(x) over the density 



Prom a different perspective, this empirical fact has recently been the object of various investigations 

6 Unless, as mentioned before, one starts at t = already with a coherent state |£, <£;iV), i.e. in a 
superfluid phase. 
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matrix p obtained from the starting state p through the action of a (completely positive) 
map; explicitly: 



p{ t) ^ p( t) = J d^) J V ^ $ ^ t) p(t) V(£, 0; N, t) , (21) 

with 

V{1 <p; N, t) = % 0; N, t){£ &N,t\; (22) 

it is a realization of a Positive Operator Valued Measure (POVM) [T5]-|17j. 

Notice that the mean value for the atom density given in f|T9|) obtained through the POVM, 
i.e. 

n p (x,t) =Tr h(x) p(t) , (23) 

is clearly distinct from that obtained using the definition ( ITTj) ; nevertheless, they both repro- 
duce the profile ffTB"]) in the case of a superfluid initial state, p = |£, <p; iV")(£, 0; N , thanks 
to the (large N) orthogonality of the coherent states. On the other hand, in the case of a 
Fock state, p^ = \k; iV)(A;; N\, the explicit evaluation of ( 1231) gives 

AT M 

which differs from the result (ITS]) obtained applying fll7p by an overall normalization factor, 
and in the weight assigned to the contribution of every single lattice site. Although in line 
of principle these differences might have experimental relevance, their actual detection looks 
technically problematic, except perhaps in the case of a double-well potential (M = 2) [19J. 
Instead, as we will see in the next Section, density correlations appear more suitable for 
studying the predictions of the POVM prescription. 



5 Correlation functions 

Besides for density estimations, absorption images can also be used to extract information 
on density correlations: by analyzing the absorption figure in distinct points one can study 
correlations in atom positions. Let us then introduce the two-point correlation function as 
the average of the following two-point operatoiQ 

h(x, x') = ijj\x)ip\x')ijj(x)ip(x') , (25) 

that, according to the two measuring interpretations discussed in the previous Section, can 
be either computed using the analog of the standard trace formula (ELZj), 

n p (x,x',t) = Tr[n(x, x') p(t)] , (26) 

7 It differs from the density-density correlation operator by a 8(x — x') contribution; on average, this term 
is suppressed by a factor 1 /N and therefore negligible in the large N limit [24] . 
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or by assuming that a single shot picture effectively projects the system state into a coherent 
one, so that operator averages should be taken using the transformed density matrix p(t) of 
( I2T1) . thus giving 

n p {x, x', t) = Tr [n{x, x') p{t)] . (27) 

In order to make an easier comparison with the actual experimental data, a further in- 
tegration with respect to the barycenter coordinate, R — (x + x')/2, is usually performed, 
so that the resulting integrated correlation function depends only on the relative coordinate 
r = x' — x [21]. Then, with a suitable normalization, one is lead to study the behaviour of 
the following two functions 

JdRn p (R-^R+^t) 
Gp[T ' t] = JdR n p {R-^t) n p (R + $,t) ' (28) 

and, by replacing the standard averages with the POVM ones, 

jdflfi,(fl-j,fl +i ,t) 

y ' { • 1 ~ /cLRfi„(fl-§,t) h„(R + l,t) ■ { ' 

These functions measure the conditional probability of finding two atoms in points separated 
by a distance r, averaged over all positions. In absence of correlations, they take a constant 
value equal to one, while values greater than one would signal the tendency of the atoms to 
aggregate, a typical behaviour for bosons. As in the case of the average density discussed in 
the previous Section, the expressions ( 1281) and ( 1291) take a particularly simple and compact 
form when the system is initially prepared in a Fock state, p = | k; N}(k; N\. Indeed, recalling 
( II Op . that gives the free evolution of the Wannier functions, one explicitly finds 



while 



5M = + WTm L TW - T) ± (h + l )(k , + De^j , (31) 



where, as before, k{ represents the initial occupation number of the i-th lattice site and 
Q = md/ht. It turns out that for a given initial state, the overall profile of (1301) and (|3~T1) 
is very similar: both functions are essentially constant and close to one almost everywhere, 
except at the origin and at positions multiple of 2ir/Q where sharp peaks occur. 

Nevertheless, the two expressions do differ in the normalization of the oscillator terms, as 
well as in their dependence on the occupation numbers. One can check that the effects of 
these differences become more and more visible as the initial configuration differs from that 
described by a balanced Fock state, the one with an equal number N/M of atoms in each 
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lattice site. Indeed, for such a state, up to terms of order l/N, one finds that both functions 
(BOD and (ED) reduce to: 



In order to appreciate the differences between the two theoretical predictions Q p (r,t) and 
Q p (r,t), one has therefore to prepare the system in an unbalanced Fock state. The way 
in which this can be experimentally realized is by superimposing a second weaker periodic 
potential to the original one; this is done through the introduction of a second laser directed 
along the optical lattice, whose wavevector k<i is chosen to be incommensurate with respect 
to Hi = ir/d, the wave vector of the first, original one [20]-[23]. When the amplitude of the 
second laser is weak, thus introducing only small perturbations to the original potential, the 
dynamics of the atoms in the lattice can still be described by the Bose-Hubbard Hamiltonian 
(CQ), where however the site depth energy e« gets a further site dependent contribution of order 
V2 sin 2 (i/? 2 7r/fi;i), proportional to the strength V2 of the second laser. 

Then, in order to prepare the system in a Fock state, one gradually increases the intensity 
of the lasers, so that the hopping term in the Hamiltonian (TjQ) becomes negligible. In this 
way, one drives the system into a ground state of the form |fc;iV), characterized by the 
filling configuration k, minimizing the total energy E = £V ejkj + ^ J^ . kj(kj — 1) with the 
constraint Ylj kj = for kj integer. The advantage of working with a bichromatic lattice is 
now apparent: it allows to change the distribution of the N atoms in the M wells by varying 
the amplitude V2 of the second potential. 

For a given value of the characteristic, physical parameters entering the Hamiltonian (TjQ), 
the configuration k that minimizes the energy can be efficiently obtained using a numerical 
simulation. We have used a Monte Carlo method implementing a simulated annealing algo- 
rithm, that better conforms to the experimentally adopted procedure, since it considers all 
atoms in the lattice at oncej^l Indeed, in a typical experiment, the atoms are first cooled in an 
harmonic potential, at the center of which the optical lattice is then slowly raised. Initially, 
the hopping dynamics allows a redistribution of the atoms in the various sites; however, this 
becomes highly suppressed at regime, leaving the system in a Fock ground state. 

To be as closer as possible to an actual experimental situation, we have chosen to work 
with the physical parameters that define the apparatus described in Ref. [2T]. Although the 
principal lattice is in that case three-dimensional, the second, weaker potential is switched 
on only along one direction, making the whole system effectively behaving as a collection of 
independent, separate one-dimensional bichromatic lattices. The whole system is filled with 

8 This simulation procedure should be contrasted with the one adopted in [52], where the N atoms are 
inserted in the lattice one by one, while minimizing the energy at each step. 




(32) 



(33) 
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about 3 x 10 5 particles, but any single one-dimensional lattice is formed by about M = 130 
sites, filled with roughly N = 170 ultracold atoms. 

The result of the simulation is summarized in Figure 1: it shows the plots of the corre- 
sponding correlation functions Q p (r,t) (green line) and G p (r,t) (red line). As expected, these 
functions are almost everywhere equal to one, except for the presence of periodic bumps. 
The difference between the two is particularly visible in the lower, secondary peaks, whose 
height is highly suppressed when the correlations are computed using the POVM prescrip- 
tion. This difference become more and more evident as the strength of laser giving rise to 
the second potential is increased, as clearly shown in Figure 2. As a result, the suppression 
of the secondary peaks can be made visible well beyond any statistical error, thus becoming 
experimentally relevant. 



G(r) with and without POVM 
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Figure 1: Behavior of Q(r, t) (green line) and of Q(r,t) (red line) with N = 170, M = 130, V2 = 
h x 9.9 kHz. 



6 Outlook 

The measuring procedure commonly used in experiments with ultracold gases, consisting 
in extracting density profiles from absorption images taken after the release of the optical 
lattice, suggests a theoretical interpretation in terms of generalized quantum measurement 
processes. The experimental evidence regarding the presence of interference fringes in single 
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height of the secondary peak 
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Figure 2: Height of the secondary peaks as a function of the intensity of the secondary laser with 
(red line) and without POVM (green line). 



shot absorption pictures irrespective from the initial state of the system suggests the use of 
a POVM based on coherent-like, fixed phase states. Within this framework, the averages 
of system observables in general differ from those obtained through mean values of the 
corresponding operator in the system state, the usually adopted paradigm in interpreting 
experimental data. 

As discussed in Section 4, these differences are hardly visible in the case of density mea- 
sures, since they are of the order of the inverse total number of atoms in the sample. Instead, 
the situation appears quite different for density-density correlations. As explicitly discussed 
in the previous Section, by preparing the system in a suitable Fock state, the profile of the 
integrated and normalized correlation function is characterized by a double series of peaks, 
the smaller of which appear to be much lowered when the correlations are computed using 
the POVM prescription. This effect is very pronounced and surely well beyond statistical 
errors for situations close to the actual experiment: this may open the way to a direct test 
of the POVM assumption. 
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